clear all
close all

t = linspace(0,2*pi,1001);
S = 1:1001;
for j = 1:1001
    S(j) = 1/2 + 2 * sommeFourier(t(j),200) / pi;
end

figure(1)
plot(t, S, 'b');